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Abstract 

We calculate the effective Polyakov line action corresponding to SU(2) lattice gauge theory on a 16 3 x 4 
lattice via the "relative weights" method. We consider a variety of lattice couplings, ranging from j3 = 1.2 in 
the strong-coupling domain, to j8 = 2.3 at the deconfinement transition, in order to study how the effective 
action evolves with j8 . Comparison of Polyakov line correlators computed in the effective theory and the 
underlying gauge theory is used to test the validity of the effective action for j3 > 1.4, while for j3 = 1.2, 1.4 
we can compare our effective action to the one obtained from a low-order strong-coupling expansion. Very 
good agreement is found at all couplings. We find that the effective action is given by a simple expression 
bilinear in the Polyakov lines. The range of the bilinear term, away from strong coupling, grows rapidly in 
lattice units as j5 increases. 
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I. INTRODUCTION 



In a recent article [1] we have applied a technique, which we call the "relative weights" method, 
to determine the effective Polyakov line action (PLA) corresponding to an SU(2) lattice gauge 
theory. The effective Polyakov line action Sp is defined as the action which results from integrating 
out all gauge and matter degrees of freedom in the lattice gauge theory, whose action is denoted 
Sl, under the constraint that Polyakov line holonomies are held fixed. In temporal gauge, where 
the timelike link variables are set to the identity matrix except on a single timeslice at, say, t = 0, 
we have 1 



exp 



Sp[U x ]\ = J DU Q (x,0)DU k D<l> ^Y[8[U x -Uo(x,0)]je SL ■ (D 



where denotes any matter fields, scalar or fermionic, coupled to the gauge field. 

In pure SU(2) lattice gauge theory, which is the case we will consider in this article, Sp can 
depend only on the trace of Polyakov line holonomies. Consider a Fourier expansion of the trace 

P x = -Tr[U x ] = a o + l Y,{ a 4 cos (4' x ) +b 4 sin (4' x )} ' ( 2) 

where the sum runs over all wavevectors q on a cubic lattice of volume L 3 , and a q = a- q , b q = —b- q 
are real-valued. The relative weights method allows us to compute the derivatives of the action 
with respect to any of the Fourier components 



dSp\ ( dSp 



da kJa k =a ' \ db k/b k =a 



(3) 



and from these derivatives we are able to deduce the action itself. In ref. [1] we found that the 
effective action is a simple bilinear expression in terms of the Polyakov lines: 

Sp = \c x J>2 _ 2c2 ^P x Q(x - y)P y , (4) 

x xy 

where 

Q(x-y) = { (V^,, I*"?! £ r ~ . (5) 
[ \x-y\>r max 

and is the lattice Laplacian. Our method also determines the constants c\ , and the range r max; . 
The validity of this expression was tested by comparing Polyakov line correlators computed by nu- 
merical simulation of the effective theory based on Sp, and of the underlying lattice gauge theory. 
We found accurate agreement, down to correlator magnitudes on the order of 10~ 5 , between the 
two sets of correlators. 

A limitation of the work in [1] is that it was carried out for a single gauge coupling /3 = 2.2 
at Nt = 4 lattice spacings in the time direction. An obvious question is how the effective action 
evolves as /3 varies from strong couplings to weaker couplings, up to the deconfinement phase 
transition. Answering this question will also allow us to check how the relative weights method 



For convenience, we adopt a sign convention for the action such that the Boltzman weight is proportional to exp [+S] . 
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performs over a range of couplings, rather than just a single coupling. At very strong couplings the 
effective action can be computed analytically, and only the nearest-neighbor term in the effective 
action is significant. At next to leading order, we have 2 



5p = jBp££flA+i 

x i=l 



]6p = 4 



(h(P)\ Nt 



(6) 



where N t is the lattice extension in the periodic time direction. In contrast, at /3 = 2.2, we found it 
necessary to include couplings among Polyakov lines separated by distances up to r max = 3 lattice 
units. We would like to study how the action, and in particular the range of the bilinear term r max , 
evolves from strong coupling through to the deconfinement transition. 

In this article we extend our previous work by computing Sp, via the relative weights method, 
for a variety of lattice couplings, starting from a strong coupling value of /3 = 1.2, and proceeding 
to the deconfinement transition at /3 = 2.3, N t = 4. All our numerical computations are carried out 
on a 16 3 x 4 lattice volume for the SU(2) lattice gauge theory, while simulations of Sp are carried 
on a three dimensional 16 3 lattice volume. Our main finding is that the effective action is well 
described by a bilinear form throughout the range of lattice couplings, and that the range of the 
kernel Q(x — y), beyond the strong-coupling limit, increases rapidly in lattice units with increasing 
J8- 



II. PROCEDURE 

The relative weights method allows us to calculate numerically the difference 

ASp = Sp[U£-S P [UZ] (7) 

of Polyakov line actions, evaluated at configurations U' x and U x respectively, which are nearby in 
the configuration space of Polyakov line holonomies. The method is based on the fact that, while 
it may be difficult to evaluate the integral in (1) directly, the ratio (or "relative weights") 

asp = exp[5 P [^]] 

" «p[W1] 

can be expressed in a form which is more amenable to numerical simulation. Let S r L ,S'[ represent 
the lattice action with timelike links Uq(x,0) fixed to U' x and U x respectively; these links are not 
integrated over. Then 

eASp = $DU k D<$> A 
fDU k D^> e s L 
_ fDU k D(j) cxp[S' L -S f []e s 'L 
fDU k D(j) e s L 

= (cxp[S' L -S'[]Y (9) 

2 For a higher-order computation, cf . [2] . 
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where (...) ;/ indicates that the VEV is to be taken in the probability measure 

(10) 



fDU k D(j) e s L 

As mentioned in the previous section, for the SU(2) gauge group the effective action Sp depends 
only on the trace P x = jTr[U x ] of Polyakov line holonomies, so we may consider two configura- 
tions in which the Fourier decomposition (2) of Polyakov lines in each configuration differ only by 
Aafc in the amplitude of a particular Fourier component. In that case we can estimate the derivative 
by a finite difference 

SS A . AS P 



da kJa k =a ^ a k 

The remaining Fourier components, which are the same for Polyakov lines in U' and U" , are 
derived from a thermalized lattice configuration. The procedure is to generate a thermalized con- 
figuration Up (jc, t) by the usual lattice Monte Carlo method, calculate the associated Polyakov line 
holonomies 

U x = U (x, l)U (x,2)-Uo(*M) (12) 

and carry out the Fourier decomposition of the corresponding Polyakov lines (2). Then pick a 
particular wavenumber k and set a# = 0. Denote the modified Polyakov lines as P x . We then 
construct two other Polyakov line configurations 

P x = acos(k-x)+fP x 

pj!= (a+Aa k )cos(k-x)+fP x (13) 

along with corresponding holonomies U X ,U X which give rise to these Polyakov lines. The constant 
/ ~ 1 — a, is chosen to ensure that the absolute values of P' x and P x are < 1 . We then compute 
(1 1) by the relative weights approach. 

For a detailed exposition of the relative weights method, including noise reduction and the 
precise way in which we choose U x , U x and /, the reader is referred to ref. [1]. 



HI. RESULTS 



We have carried out simulations of pure SU(2) lattice gauge theory on a 16 3 x 4 lattice volume 
at the following /3 values: 

/3 = 1.2, 1.4,1.6, 1.8,2.0,2.1,2.2,2.25,2.3 (14) 

All of these values lie inside the confined phase with the exception of j8 = 2.3, which lies essen- 
tially right at the deconfinement transition for N t = 4 (the precise transition point is /3 C = 2.2986(6) 
[3]). We also calculate the derivative (11) at a values 

a = 0.05,0.10,0.15,0.20 (15) 
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and for a range of lattice momenta < ki < 2.9, where 



\ 



4fsin 2 (k-) (16) 



i=l 



The components ki of the wavevector k are k\ = ^-nii, with L = 16 the spatial extension of the 
lattice, and m, = integers (including zero). 

The striking fact is that, in all cases, the derivative (11) of Sp is linear in a. This is evident 
when we plot 



1 1 fdSp 



aL 3 \da k 



vs. k L (17) 

a^—a 



as shown in Fig. 1 (/3 = 1.2—1.8) and Fig. 2 (/3 = 2.0 — 2.3). In these figures the data points 
displayed at kt = are actually equal to the data values divided by a factor of two, for reasons we 
will explain. The important feature to notice in each of these plots is that at any given /3 and k L , the 
data points at different a essentially coincide. This means that the first derivative of Sp is linear 
in a, and it follows that Sp itself is quadratic in a k for any k. As a consequence, Sp will also be 
quadratic in the position space variables (i.e. the Polyakov lines P x ) and we can write this effective 
action in the form (4). The problem is then to determine c\ , c% and Q(x — y) from the data. 

Let Q{k£) be the finite Fourier transform of Q(x). This leads to derivatives 

1 J dS P [U x (a k )] \ = (oc(±c l -2c 2 Q(k L )) k L ^0 ^ 
L 3 V da k J ak=a \ 2a {\ Cl -2c 2 Q{Q)) k L = 

The relative factor of two in the ki = and k^ > cases is due to the fact that 1 = L 3 , while 
£ x cos 2 (^- x) = jL 3 . The ki > data should extrapolate, as ki — > 0, to a value which is half the 
result at ki = 0. For this reason, in Fig. 1 and subsequent figures, the point shown at ki = is the 
data value at ki = divided by two. 

Let us begin with the data for lattice couplings /3 in the intermediate/weak-coupling regime 
2.0 < /3 < 2.3 shown in Fig. 2. In addition to the a-independence, what is striking about this data 
is that it is clearly linear for most of the range of lattice momenta. If the data were linear for the 
entire range, then we would have 

5/>=k£P 2 -2c 2 £P l ( V / ^Vi) P yi (19) 

x xy x y 



In this case the kernel is Q(x — y) = ( y — ^jjxy which translates in momentum space to Q(k) = k^. 
This corresponds to the linear behavior seen in Fig. 2. However, the action (19) has a long- 
range coupling between Polyakov lines which, in the first place, is inconsistent with the non-linear 
behavior seen at low ki, and in the second place would violate one of the assumptions of the 
Svetitsky-Yaffe conjecture [4], which postulates only finite-range couplings in the effective action. 
A simple finite range ansatz for Q(x — y), having the linear behavior in momentum space seen at 
high ki, is the expression (5) proposed in ref. [1]. Since Q(k) — > ki at high k^, the constants 
ci,C2 are determined from a linear fit of the data in the linear (higher- momentum) regime to the 
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FIG. 1. Derivatives of Sp with respect to Fourier components a^, evaluated at = a, and divided by a, for 
intermediate/strong-coupling values of j3 = 1.2 — 1.8. The data, obtained by the relative weights method, is 
plotted vs. kL. The point displayed at ki = is the data value divided by two, for reasons explained in the 
text. 



expression 



1 l_ fdS P [U x (a k )} 



a 



da k 



a/ l =a 



2 Cl 



2c 2 h 



(20) 



The action (4) is then well-defined, given the finite range cutoff r max . 

To determine r max , we look for the value which satisfies the lower {ki = 0) identity in (18), i.e. 



1 1 (dS P 
ccL 3 \dao 



ci-4c 2 0(O) 



(21) 



as accurately as possible. We cannot satisfy this relationship exactly because, while /3 can be 
varied continuously, r max cannot, since on the lattice it is the square root of a sum of squared 
integers. Nevetheless, if r max is not too small, we can come close to satisfying the equality. This 
is essentially an optimization problem. For each trial r max , compute Q{x - y ) from (5), and then 



3 The precise choice of interval in ki to carry out the linear fit is a potential source of systematic error in c\ ,C2- In 
practice we choose a low-momentum cutoff which minimizes the % 2 value of the linear fit. 
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FIG. 2. Same as Fig. 1, for intermediate/weak-coupling values j8 = 2.0 — 2.3. 



Fourier transform to obtain Q{ki). The optimal choice of r max is the value which comes closest to 
satisfying (21). We can then compute the right hand side of (18) at all ki, and see how well it fits 
the data. The fractional deviation between the left and right hand sides of (21) 



^(S)_- c >- 4c ^°) 



oo=a (22) 



1 1 f dS P \ 
az7 \dao J 



ao=a 
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FIG. 3. Comparison of c\ /2 — 2ciQ{k) to the relative weights data, where Q(k) is the Fourier transform of 
the kernel (5), for couplings jS = 2.0 — 2.3. Red squares are data points, blue dots are the ci /2 — IcjQik) 
values, and the green line is the linear fit used to determine c\,C2- 



is on the order of 1% or less, for /3 > 2.0, and about 2% at /3 = 2.0. 

The results, shown in Fig. 3, seem to be quite satisfactory. The green line is the linear fit used 
to compute c\,C2, the blue dots correspond to jCi — 2c2<2(k). The red open squares are the data 
points already seen in Fig. 2, this time with no distinction on a values. 

The essential test of the proposed effective action (4) with kernel (5) is the comparison in the 
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SU(2) 16 3 x4 beta=2.20 su(2) 16 3 x4 b eta=2.25 




R 

(c)jS=2.3 

FIG. 4. Comparison of Polyakov line correlators G(R) computed by simulation of the effective Polyakov 
line action, and by simulation of the underlying lattice gauge theory. 

Polyakov line correlators 

G{x-y) = (P x P y ) (23) 

computed in both the effective theory, and in the underlying lattice gauge theory. Results for 
on-axis separations at /3 = 2.2,2.25,2.3 are shown in Fig. 4. The correlators in the lattice gauge 
theory have been calculated using Liischer-Weisz noise reduction [5]. 

The agreement between the correlators in the effective theory and the lattice gauge theory is 
remarkably accurate, down to values of G(x — y) ~ 10~ 5 . However, as /3 is reduced, we find that 
the data in the effective theory seems to flatten out at around 10~ 5 , as can be seen in Fig. 5. This 
flattening may simply be a finite volume effect in the numerical simulation of the effective theory. 
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FIG. 5. Same as Fig. 4, for j3 = 2.0,2.1. Note the plateau, in the correlator of the effective theory, around 
G(R) = icr 5 . 



The argument is as follows: For separations R = \x — y\ which are several times greater than a 
correlation length, say for R > R ma x, Polyakov lines are almost uncorrelated. Let us denote the 
average value of a Polyakov line in a given thermalized configuration as 



1 



and also define 



y>R>R„ 



Ljc PxPy 0(R— R mi 
HxT.y d (R — Rmax) 



(24) 



(25) 



where (x) is the Heaviside theta function. Now for R > R ma x Polyakov lines P x and P y are essen- 
tially uncorrelated, and if the lattice volume L 3 is much greater than R max , we may approximate P x 
and P y in the double sum by their average values in the configuration. But this means that in each 
configuration 



( p x p y)R>R„ uu 



Px 



(26) 



On any given site, the typical magnitude of P x is perhaps on the order of 0.2, which means that the 
average value on the lattice, in a typical configuration, will be on the order of P x ~ ±0.2/L 3//2 . For 
L = 16, this gives us an estimate of 



y}R>R„ 



10 



(27) 



in any thermalized configuration in the effective theory on a 16 3 lattice volume. This may explain 
why the measured values of G(x — y) seem to plateau at around 10~ 5 at these lattice volumes. 

Reducing the lattice coupling below /3 = 2.0, we enter the regime of strong couplings. Our 
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FIG. 6. Comparison of the best fit c\/2 — 2c%k^ to the relative weights data at strong-couplings j8 = 1 .2, 1 .4. 



results for 

1 1 fdSp 



a L 3 \ dcik 



vs. k L (28) 

a k =a 



at /3 < 2.0 were seen in Fig. 1. In these plots there is much more evidence of curvature, and the 
part of the graph which fits a straight line disappears as /3 is reduced. It turns out that the curved 
section of the plots are compatible with Q(k) = k\, and in fact this fits the data at /3 = 1.2, 1.4 
perfectly, as we see in Fig. 6. At /3 = 1 .6, 1 .8, however, there is still a portion of the data which fits 
a straight line, and so we need an ansatz for Q(k) which interpolates between the quadratic form 
at small kj_ and linear at large k^. A simple choice is 



Q(k) = ^k 2 + m 2 (29) 

and so we do a best fit of the data at all ki to the form 

1 1 f dS P [U x (a k )} \ 1 ,f^~~ 2 „m 

™73 1 =-ci-2c 2 yk +m 2 (30) 

aU \ da k J ak=a 2 v 

The resulting fits, for lattice couplings j8 = 1.6, 1.8, are shown in Fig. 7. 

With the constants ci,C2,m 2 determined from a best fit to the data, the position space ker- 
nel Q(x — y) is determined from a fast Fourier transform of Q(k), and we can again compare 
the Polyakov line correlators computed in the effective theory with the corresponding correlators 
computed in the underlying lattice gauge theory. The results are shown in Fig. 8, this time includ- 
ing off-axis separations. Once again, the agreement is very good, although the falloff in position 
space is so rapid that the 10~ 5 plateau sets in at rather small values of |x — y\. In these figures 
the Polyakov correlators in the lattice gauge theory are computed in the standard way, without 
the Liischer-Weisz noise reduction, and these also seem to show a plateau at the larger lattice 
separations. 

Since Q(x — y) corresponding to a square root ansatz does not have a sharp cutoff in the range 
of R = \x—y\, the corresponding effective action becomes challenging to simulate numerically as 
m is reduced. We have not yet investigated this form of the action in the weak coupling regime. 
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FIG. 7. Comparison of the best fit (blue line) of the square-root ansatz c\/2 — 2c2 y k[ + m 2 to the relative 

weights data at strong/intermediate-couplings j8 = 1.6,1.8. Also shown is a fit to ci/2 — 2c2&£ at lower 
momenta (green line). 
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FIG. 8. Comparison of Polyakov line correlators G(R) computed in the effective theory (green dots) and in 
the underlying lattice gauge theory (red squares), at the intermediate/strong couplings /3 = 1.6, 1.8. These 
plots include off-axis separations. 



Finally, at our two strongest couplings, /3 = 1.2 and 1.4, there is no need for an interpolating 
form, since the Q(k) = k\ ansatz fits the data quite well, as seen in Fig. 6. The data is fit to the 
form 

1 1 ( dSp[U x (a k )} \ 1 2 

M^^U = 2 C '- 2CA <31) 
and a short calculation gives us the effective action 

3 i 

Sp = 4c 2 ££/yWi+ (-ci - 12c 2 )£p2 (32) 

x i—1 ^ x 
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J_ 

1.2 
1.4 
1.6 

1.8 

2.0 
2.1 
2.2 
2.25 
2.3 



G(*-y) 



xy 



xy 



-V? + m 2 



eq. (5) 
eq. (5) 
eq. (5) 
eq. (5) 
eq. (5) 



xy 



xy 



0.1714(2) 
0.3196(5) 
4.10(7) 

1.969(8) 

2.93(1) 
3.63(1) 
4.417(4) 
4.70(1) 
4.812 a 



C2 

0.007148(8) 
0.01327(2) 
0.219(2) 

0.269(1) 

0.313(2) 
0.397(2) 
0.498(1) 
0.541(2) 
0.563(2) 



m 



4.01(4) 
2.37 (2) 



V3 
V5 
3 

To 
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a The actual constant derived from the fit to the linear portion of the data was 4.77(1). However, the correlator at 
j3 = 2.3 is extremely sensitive to value of c\, and a small adjustment to 4.812 greatly improves the agreement with 
the Polyakov line correlator derived from lattice gauge theory. 

TABLE I. Constants denning the effective Polyakov line action (4) for pure SU(2) lattice gauge theory on a 
16 3 x 4 lattice. 



The best fits to (31) give the constants 

ci =0.3196(5) c 2 = 0.01327(2) = 1.4 



c 2 = 0.1714(2) c 2 = 0.007148(8) j8 = 1.2 (33) 

With these numbers, we find that 

|ci -12c 2 = (-0.56 ±0.34) x 10~ 3 (j8 = 1.4) 

| C1 -12c 2 = (-0.76 ± 1.4) xlO- 4 (j8 = 1.2) { ' 

which are essentially consistent with zero. Then, comparing the relative weights and strong- 
coupling results for the PLA at /3 = 1.4, we find 

0.05308(8) P x P x +i relative weights 

S P ={ Q3 = 1.4) (35) 



0.0522 E?= i PxPx+t strong coupling 



while at j3 = 1.2, 



0.02859(3) ZxZl=iPxP x +i relative weights 
S P ={ OB = 1-2) (36) 

0.02850 LxL]=i PxPx+i strong coupling 

In both cases there appears to be good agreement between the strong-coupling expansion and the 
relative weights result. 

To summarize, the effective Polyakov line action is given by (4) with the parameters and kernel 
Q(x — y) listed in Table I. Note that the range of the kernel in lattice units rises from r max = \/3 ~ 
1.73, at /3 = 2.0, to r max = \/13 w 3.61 at the deconfinement transition. 
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IV. CONCLUSIONS 



We have calculated the effective Polyakov line actions corresponding to pure SU(2) lattice 
gauge theories on a 16 3 x 4 lattice volumes in an interval of lattice couplings ranging from /3 = 
1.2, which is deep in the strong-coupling regime, up to /3 = 2.3, which is at the deconfinement 
transition. At each coupling, the effective lattice action has a simple bilinear form (4), although the 
range of the bilinear kernel Q(x — y) varies from only nearest neighbor couplings, in the strong- 
coupling regime, up to separations of |jc — y \ = 3.61 lattice units at the deconfinement transition. 
This extends the work of ref. [1], where only the lattice coupling /3 = 2.2 was considered. Our test 
for checking the validity of these effective actions is the comparison of Polyakov line correlators 
calculated in the effective theory and the corresponding lattice gauge theory. In every case, the 
comparison works out quite well, down to correlator values on the order of 10~ 5 . There have 
been other approaches to calculating the effective Polyakov line action, including strong-coupling 
expansions [2], the Inverse Monte Carlo method [6], and the Demon approach [7], resulting in 
effective actions of varying complexity, but we do not believe that these have yet demonstrated a 
comparable agreement in Polyakov line correlators at the larger /3 values, at least not beyond two 
or three lattice spacings. 

So far we have not tried to calculate the effective action inside the deconfmed regime. Pre- 
liminary work suggests that the range of the kernel in this region region is rather large in lattice 
units, which makes the effective action difficult to simulate, and also it may be necessary to intro- 
duce higher powers of P x in the potential of the effective theory. We leave this problem for future 
investigation. 

Although we believe that deriving the effective Polyakov line action is of interest in itself, our 
ultimate goal is to apply our approach to the sign problem [8] . There is no sign problem in SU(2), 
nor is there a sign problem in any pure gauge theory, so the next step in our program will be to 
apply the relative weights method to the SU(3) gauge group, first without and then with matter 
fields. Our method only supplies the effective action at zero chemical potential. However, as 
explained in ref. [1], the effective action corresponding to a lattice gauge theory at finite chemical 
potential may be obtained from the effective action at zero chemical potential by a simple change 
of variables. If our method is successful for the SU(3) gauge group, then the strategy would be 
to apply one or more of the methods in refs. [2, 9-11], which were developed to tackle the sign 
problem in Polyakov line actions, to the problem of determining the phase structure of the theory. 
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